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ABSTRACT 

We investigate whether the large scale structure environment of galaxy clusters imprints a 
selection bias on Sunyaev Zel’dovich (SZ) catalogs. Such a selection effect might be caused 
by line of sight (LoS) structures that add to the SZ signal or contain point sources that disturb 
the signal extraction in the SZ survey. We use the Planck PSZl union catalog ( [Planck Collab-| 
oration et al.|2013^ in the SDSS region as our sample of SZ selected clusters. We calculate 
the angular two-point correlation function (2pcf) for physically coii'elated, foreground and 
background structure in the RedMaPPer SDSS DR8 catalog with respect to each cluster. We 
compare our results with an optically selected comparison cluster sample and with theoretical 
predictions. In contrast to the hypothesis of no environment-based selection, we find a mean 
2pcf for background stnictures of -0.049 on scales of < 40', significantly non-zero at ~4cr, 
which means that Planck clusters are more likely to be detected in regions of low background 
density. We hypothesize this effect arises either from background estimation in the SZ survey 
or from radio sources in the background. We estimate the defect in SZ signal caused by this 
effect to be negligibly small, of the order of ~ 10 of the signal of a typical Planck detection. 
Analogously, there are no implications on X-ray mass measurements. However, the environ¬ 
mental dependence has important consequences for weak lensing follow up of Planck galaxy 
clusters; we predict that projection effects account for half of the mass contained within a 15’ 
radius of Planck galaxy clusters. We did not detect a background underdensity of CMASS 
LRGs, which also leaves a spatially varying redshift dependence of the Planck SZ selection 
function as a possible cause for our findings. 

Key words: galaxies: clusters: general - cosmology; observations. 


1 INTRODUCTION 

Clusters of galaxies play a major role in astrophysics and cosmol¬ 
ogy, as they can be used to put constraints on the dark matter con¬ 
tent of the universe. Furthermore galaxy clusters are particularly 
sensitive to the interplay of dark matter and dark energy. They are 
cosmological probes that could potentially help to distinguish be¬ 
tween dark energy and modified gravity explanations for the ac¬ 
celerating expansion of the universe (for a review, see [Allen et ah] 
[2011[[Borgani & Kravtsov|2011 [[Weinberg et al.|2013^. 

A variety of different methods for cluster detection and mass 
measurement exists. Gravitational lensing probes the dark and lu¬ 
minous matter distribution of a cluster by measuring the distortion 
of background galaxies (weak lensing, for example in [Hoekstra] 
[et al.|200T][Gruen et al.|2013[|2014^ , or by detecting multiple im¬ 
ages of single background galaxies close to the LoS of the cluster 
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core (strong lensing, for example in [Zitrin et al.|2012j[Eichner et al.[ 
[2013[[Monna et al. [2014[ l. The most widespread method for optical 
cluster detection, the so-called red sequence method jGladders &[ 
[Yee|2005[[Koester et al.|2007[[Rykotr et al.|2014^ is based on spatial 
overdensities of red galaxies. Further methods include the observa¬ 
tion of the X-ray Bremsstrahlung emission by the hot gas in the 
intra-cluster medium (ICM, [Piffaretti et al.||20lT[ [Vikhlinin et al.| 
[2009[ [Mantz et al.|2010) and the observation of inverse Compton 
scattering of the cosmic microwave background (CMB) photons by 
the ICM, which is known as the SZ effect jSunyaev & Zeldovich[ 
[1972[ >. The latter describes the distortion of the CMB spectrum 
along the LoS through clusters and groups. The amplitude of the 
SZ effect is proportional to the dimensionless Compton parameter 
y, defined as the integral over the thermal electron pressure along 
the LoS: 

y=-^ f Pdl, (1) 

while the integral over a solid angle yields the SZ observable Y: 
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DlY = Dl fydQ=-^ f PdV, (2) 

J J 

where ctt is the Thomson cross section, the rest energy 
of the electrons and the angular diameter distance. 

All of these methods may have selection effects induced by 
structures along the LoS. Lensing, for example can yield biased 
mass estimates when there are groups along the LoS which con¬ 
tribute to the shear signal (e.g. |Spinelli et al.|2012) . The X-ray sig¬ 
nal of LoS structures can stack, resulting in a biased mass estimate. 
The same is true for the SZ effect, however more severely as the 
SZ signal is proportional to the gas density p, while the X-ray flux 
is proportional to p^, making the effect of LoS structure on SZ sig¬ 
nals much larger at larger angular separation. Due to this reason we 
will investigate whether SZ selected clusters are possibly biased 
by structures along the LoS, either by physically uncorrelated fore¬ 
ground or background structures, or by correlated structures at the 
same redshift as the cluster itself. 

Several effects could potentially contribute to a selection bias. 
The blending of the SZ signal of the detected cluster with groups 
along the LoS could bias the SZ estimate high and cause clusters 
along overdense lines of sight to be more likely detected. If, on the 
contrary, unresolved groups in the vicinity of clusters increase the 
background level, this could lead to a lower detection probability 
as the signal from the cluster is partly suppressed by the wrong 
background estimate. Furthermore, if the background of a cluster 
is contaminated with radio-loud galaxies, this could raise the noise 
such that clusters with a weak SZ signal are not detected. 

In this paper we address this question by analyzing the pro¬ 
jected group environment of SZ-selected clusters from the Planck 
PSZl union catalog (Planck Collaboration et al.||2013a[ l and test 
for group overdensities or underdensities along the LoS in the fore¬ 
ground, background and at the redshift of the clusters. The group 
sample is taken from the RedMaPPer red-sequence catalog based 
on SDSS DR8 photometry (Rykoff et al.||2014| |Rozo & Rykoff| 
|20T4| [Rozo~erar]|20T4) . We compute the angular two-point cor¬ 
relation function (2pcf) of galaxy clusters and groups for differ¬ 
ent subsamples of our catalogs (correlated, foreground and back¬ 
ground structures) to quantify correlated and physically uncorre¬ 
lated group overdensities and underdensities. We compare these 
results to the 2pcf obtained for an independent cluster sample with 
similar redshift and richness distribution, drawn as a subsample of 
the RedMaPPer SDSS DR8 catalog, and to theoretically predicted 
values. 


1.1 Motivation 

We briefly discuss several possible effects that could cause a se¬ 
lection bias. The filter function that is used for the Planck cluster 
detection might estimate a too large background value if there are 
groups surrounding the cluster that contribute to the signal, which 
could lead to a decreased detection probability in crowded fields 
as the subtracted background estimate is too large. On the other 
hand, the clusters are detected by combining six frequency bands 
with different filter sizes, so it is rather unlikely that this still causes 
problems when detecting clusters based on the differential signal. 

Another possible origin of a selection effect might be radio- 
loud galaxies in the background. |Donoso et al.| ( |2010[ > state that 
radio-loud active galactic nuclei (RLAGNs) are predominantly 
found in dense environments compared to radio quiet galaxies 
and regular red luminous galaxies (LRGs) at redshifts 0.4 < z < 
0.8. They conclude that this clustering effect is stronger for more 


massive RLAGNs. In [Yates et al.| ( [T989[ l the clustering effect of 
RLAGNs at z 0.2 is compared to the one at z ~ 0.5, with the 
result that the latter objects are found in environments three times 
denser on average. They also state that more powerful RLAGNs are 
found in denser environments than less powerful ones. Based on 
these findings we hypothesize that a high background group den¬ 
sity entails a higher probability of containing radio sources and thus 
increases the noise along the line of sight, potentially leading to a 
lower SZ detection probability for clusters in dense background en¬ 
vironments. 

This paper is structured as follows. In section we describe 
the Planck PSZl union catalog and the RedMaPPer SDDS DR8 
group catalog, as well as our matching of these two. In section]^ 
we briefly discuss fwo-point correlation functions. Furthermore we 
describe our method of generating random points for the Planck 
catalog and the procedure of defining the cluster comparison sam¬ 
ple out of the RedMaPPer catalog. We also include the description 
of our theoretical prediction of the 2pcf. In section]^ we present 
our results, give a detailed description of our error estimation and 
generalized x~ analysis and we estimate the implications of the 
measured effect on SZ and lensing analyses of Planck clusters. We 
conclude in section|3 


2 DATA 

2.1 The Planck PSZl catalog 


The Planck PSZl union catalog is a cluster catalog, covering the 
whole sky based on SZ detections using the first 15.5 months of 
Planck survey observations. It contains a total of 1227 clusters, 861 
of which are confirmed while fhe remaining 366 are cluster candi¬ 
dates (Planck Collaboration et al.|2013a| >. The Planck satellite fea¬ 
tures a low frequency and a high frequency instrument, the former 
covers the bands at 30, 44 and 70 GHz (Planck Collaboration et al.j 
| 2013ej while the latter operates at frequencies of 100, 143, 217, 
353, 545 and 857 GHz (Planck Collaboration et al.||2013b| > with 
angular resolutions between 9.53’ and 4.42’ FWHM, for a total 
of nine detection bands. The channel maps of the six highest fre¬ 
quency bands (100 to 857 GHz) were used to build the SZ-detection 
catalog, in order to avoid problems caused by strong radio point 
sources in cluster centers, which typically have steep spectra and 
thus do not appear in the high frequency bands (Planck Collabora^ 
|tion et al.|2013a[ l. 

The generalized NFW (Navarro et al.|19 ^ profile from |Ar-| 
|naud ef aT](2010^ was adopted for fhe cluster detection. 

Three detection algorithms were used to create the cluster cat¬ 
alog, two realizations of the Matched Multi-filter (MMF) method 


(|Herranz et al. 

2002 

Melin et al.|2006|l 

ICarvalho et al. 

2009 

201211. 


The MMF method detects clusters by using a linear com¬ 
bination of maps and a spatial filtering to suppress foregrounds 
and noise. The two implementations (MMFl and MMF3) split the 
whole sky in 640 patches of size 14.66 x 14.66 square degrees cov¬ 
ering 3.33 times the area of the sky {MMFl), and in 504 patches of 
size 10 X 10 square degrees covering 1.22 times the area of the sky 
{MMF3). The MMF3 algorithm is run in two iterations: the second 
is centered on the positions of the candidates from the first one, 
rejecting all candidates that fall below the signal-to-noise (S/N) 
threshold. The matched multi-frequency filter optimally combines 
the six frequencies of each patch and the resulting sub-catalogs for 
all patches are finally merged together to a single SZ-catalog per 


© 0000 RAS, MNRAS 000. 000-000 









































Environment-based selection effects of Planck clusters 3 


method, selecting the candidate with the highest S/N ratio. For es¬ 
timating the candidate size, the patches are filtered over the range 
of potential scales, selecting the scale with the highest S/N of the 
current candidate. Finally, the SZ-signal is estimated by mnning 
MMF with fixed cluster size and position. 

Powell Snakes is a Bayesian multi-frequency detection algo¬ 
rithm, optimized to find compact objects in a dilfuse background. 
After cluster detection, PwS merges all intermediate sub-catalogs. 
The cross-channel covariance matrix is calculated directly from the 
pixel data, which is done in an iterative way to minimize the con¬ 
tamination of the background by the SZ signal itself. In each iter¬ 
ation step, all detections in the same patch with higher S/N than 
the current target are subtracted from the data before re-estimating 
the covariance matrix. This so-called “native” mode of background 
subtraction produces S/N values 20% higher than those of the MMF 
method. In order to emulate the estimation of the background noise 
cross-power spectrum of the MMF method, PwS is run in “compat¬ 
ibility” mode, skipping the re-estimation step. 

Each of the three detection algorithms creates a catalog of SZ 
sources with an S/N ratio ^ 4.5. Obvious false detections are re¬ 
moved from each of the three individual catalogs ( [Planck Collabo- | 
[ration et al.|2013^ . 

The union catalog contains all sources that have been detected 
by at least two algorithms with S/N > 4.5 within a distance of 5', 
fixing the position of the MMF3 detection or, in case of no MMF3 
detection, keeping the position of the PwS detection. 


2.2 The RedMaPPer SDSS DR8 catalog 

The Red Sequence Matched-filter Probabilistic Percolation 
(RedMaPPer) algorithm ([R ykoff et al.[2014) is a red-sequence clus¬ 
ter finder based on the optimized richness estimator A jRykoff et al.[ 
\2Q\2\ . It has excellent photo-z performance and A has been de¬ 
signed to be a low-scatter mass proxy ([Rozo & Rykoff|2014||Rozo| 
[et al.|20i4) . The algorithm is divided into two stages. The first is a 
calibration stage where the red sequence model is derived directly 
from the data by relying on spectroscopic galaxies in galaxy clus¬ 
ters: given an initial model of the red-sequence, RedMaPPer selects 
cluster member galaxies, uses these to derive a new red-sequence 
model, and then iterates the whole procedure until convergence is 
achieved, as which point the red-sequence model is adequately cal¬ 
ibrated. The second is the cluster-finding stage, where RedMaPPer 
utilizes the red-sequence model to search for clusters around every 
galaxy in the SDSS. This work uses the updated version (v5.10) of 
the original RedMaPPer catalog of[Ry koff et al.[ ( [2014^ presented 
in [Rozo et al.[j2014l l. 


2.3 Matching of the Planck and RedMaPPer catalogs 

In order to match the Planck PSZl union catalog to the RedMaPPer 
SDSS DR8 cluster catalog we used an algorithm similar to the one 
described in [Rozo et al.[ ( [2014T l. We find all matches in the RedMaP¬ 
Per catalog in a radius of 10' around each Planck cluster. In the 
case of multiple matches we define the best match as the RedMaP¬ 
Per system with the highest richness. We flag all matches with a 
redshift difference between the Planck and RedMaPPer redshift of 
more than 3cr, where cr corresponds to the redshift error given in 
the RedMaPPer catalog. This gives us a total of 290 matched clus¬ 
ters. 



Figure 2. Log-log histogram of the detection probability. 


Outlier rejection 

We reject all matches that are obvious SZ-projections (5 cases), as 
identified by [Rozo et al.[j20l4| l. All clusters that have been flagged 
as 3cr redshift outliers are cross-matched with the [Rozo et al.[p0l4| l 
table of redshift outliers, and in the case of an incorrect Planck red¬ 
shift and a correct RedMaPPer redshift, we accept the cluster using 
the RedMaPPer redshift and vice versa. Furthermore we reject clus¬ 
ters with bad z-matching when a visual inspection identified them 
clearly as a mismatch (one case only), and we reject clusters that 
are outliers in the mass- fsz -plane (according to [Rozo et al.|2014l l 
due to a low RedMaPPer richness (one case only). After rejecting 
all outliers the final matched catalog includes 265 clusters. 


3 METHODS 

3.1 Two-point correlation function 

We measure the crowding of clusters and groups with the angular 
two-point correlation function, which traces the amplitude of clus¬ 
ter/group clustering as a function of their separation. The angular 
correlation function w(9) is defined as the excess probability over a 
random, uncorrelated distribution of finding two objects separated 
by an angle 6 . The probability of finding two objects in two in¬ 
finitesimal solid angle elements SQi and SQ 2 separated by angle 0 
then reads: 


dP = nin 2 (1 -I- w(9)) 6 Q. 16 Q. 2 , 


(3) 


with III and ^2 being the mean cluster/group densities in both sam¬ 
ples. 

A multitude of different estimators exist for calculating the two- 
point correlation function from data catalogs. [Kerscher et al.[ 
( [2000^ , who compared the nine most important of these estima¬ 
tors in terms of the cumulative probability of returning a value 
within a certain tolerance of the real correlation, show that the 
Landy-Szalay estimator (hereafter LS, [Landy & Szalay|1993^ per¬ 
forms best according to their criteria. Hence we adopt this estima¬ 
tor, which reads: 


itiLsiS) = 


DD - 2DR + RR 
~RR 


(4) 
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Figure 1. Redshift (left) and richness distribution (right) of the Planck sample (black) and the comparison sample (red). Planck error bars are Poissonian. 
Comparison sample error bars are not shown but are of comparable size. 


or, in case of a cross-correlation between two different sam¬ 
ples: 

„ DiDi —D^Ri —D2R1 + R1R2 

uiLsiS) = --> (5) 

K1K2 

where DD, DR and RR stand for the data-data, data-random 
and random-random pair counts, respectively. All pair counts in eq. 
[^are normalized to the total number of data pairs in the respective 
samples. Random points account for geometrical effects like sur¬ 
vey boundaries and masks in the survey area. We do not want the 
random points to correct for environment-based detection effects, 
since this is the effect we want to measure, so we are using random 
points where the true detections have been erased. The pair counts 
have been computed using the 2d-tree code Athena (|Kilbinger et al.| 
|2m4l l. 

3.2 Generation of random points for the Planck catalog 

The LS estimator (eq.|^ needs a random catalog for each data cat¬ 
alog, in order to correct for geometrical effects that could mimic a 
signal. 

There are two effects that might imprint a spatial variation on 
the Planck detection function: the variation in the noise level and 
the distance from the galactic disk. In this section we describe our 
approach to generate random points for the Planck catalog taking 
into account the varying noise level. The variation of the detection 
probability as a function of distance from the galactic disk is inves¬ 
tigated in appendix [a] 

Since the noise level of the Planck observations varies over 
the SDSS region, we need to test whether the density of SZ detec¬ 
tions has a significant correlation with the noise level that has to 
be accounted for when generating a random catalog. We use the 
Planck SMICA map (which comes in Healpix jGorski et al.|2005) 
coordinates with =2048, which is 50331648 pixels, resolution 
~ 1.7'), which uses an optimal combination of the nine frequency 


bands ( [Planck Collaboration et aH2013c^ to generate a map dis¬ 
playing the weighted average noise of all channels, averaged to 
3072 pixels, to find the noise at the position of each cluster. 

We test for correlation of the density of Planck detections 
with the noise quantitatively. In this case we assume the number 
of Planck detections per unit area / to be a power law of the noise 
N per redshift bin i with redshift dependent exponent a, : 

/■(A) = /;(!)-A"'. (6) 

We perform a likelihood analysis over the parameters /i( 1) and 
Ui, by calculating the expected number of clusters in each sky cell 
via eq. and computing the Poisson probability with the actual 
number of detections in that sky cell. The power o-, scatters around 
and is consistent with zero for redshifts z < 0.5. Above this redshift, 
we find a, x: 0.8. In conclusion, the noise level has no impact on 
detections for z < 0.5. We decide to remove all clusters with z > 0.5 
from our catalog, bringing our sample size down to 250 clusters. 
Based on these findings, we decide to use uniformly distributed 
random points for the Planck catalog. 

We use the Planck survey mask (Healpix =2048) to de¬ 
fine the region where to generate the points and cut them afterwards 
to the SDSS footprint. The random points are generated in Healpix 
coordinates to ensure a uniform distribution over the sky. 

3.3 Generation of random points for the RedMaPPer catalog 

To generate a random point catalog for RedMaPPer, we first draw 
a random position in the sky, and then randomly draw a RedMaP¬ 
Per cluster. Given the assigned cluster redshift and richness, we 
use our cluster model to randomly draw cluster galaxies to create 
a synthetic cluster. We then run RedMaPPer at this location, and 
determine whether the synthetic cluster is detected of not. The pro¬ 
cedure is repeated 100 times, and we calculate the fraction of times 
w that the cluster was detected at that location. The quantity w is 
the weight assigned to this random point. 
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Figure 3. Sketch of the selection method for the group catalogs. The green, 
blue and red volumes show the selection for the correlated, foreground and 
background samples respectively. The total z-depth of the green volume is 
0.1, with the cluster in the center. The angular radius for all volumes is 1 ° 
(see cyan circles). 


There is one subtlety associated with the above procedure: by 
random luck, some fraction of our synthetic clusters will overlap 
with real RedMaPPer clusters in both location in the sky and red- 
shift. If one did not remove the galaxies associated with the origi¬ 
nal RedMaPPer cluster before placing the synthetic cluster at that 
location, upon running RedMaPPer one will always find a cluster 
there (i.e. the original cluster), and one would erroneously conclude 
u; = 1 irrespective of the details of the synthetic cluster. Thus, it is 
critically important to remove the original RedMaPPer galaxy clus¬ 
ters from the galaxy catalog prior to drawing our random points. We 
erase clusters probabilistically: given a cluster of richness lambda 
at redshift z, we collect all of its member galaxies, and remove each 
galaxy according to the assigned membership probability, so that a 
galaxy that is 90% likely to be a cluster member is removed from 
the galaxy catalog with 90% probability. 


3.4 Definition of a comparison sampie 

We want to test whether SZ selected clusters are generally found in 
a different environment than similar (in terms of redshift and rich¬ 
ness) clusters that are selected for their optical properties. Thus we 
need to compare the cluster-group two-point correlation functions 


that we obtain for groups in the vicinity of Planck selected clusters 
to an independent sample of optically selected clusters that resem¬ 
bles the selection function of our main sample in terms of their red¬ 
shift and richness distribution. To this end, we need to model the 
Planck detection probability. We assume that the probability that a 
RedMaPPer cluster is detected by Planck takes the form: 


f^det — 


1 

2 


1 -t- erf 


R ^det \ 

dla ] 


(7) 


where erf is the error function, Ajet is the richness at which the 
detection probability is 50% and cr the scatter in richness at fixed 
SZ signal. Pdet states the probability that a cluster of given richness 
A is detected by the Planck survey, if it was inside the survey area. 
We use the RedMaPPer SDDS DR8 catalog, calculate P^a for each 
cluster and assign it as a weight to the cluster itself. 

We parameterize the redshift evolution of Ajet and a as: 


Adet = a 


( 8 ) 


and 


cr = aM +zff. 


(9) 


To find the optimum values for and ft a-, we perform 

a likelihood analysis in these four parameters: 

ln(L) = In [Pdet(0] + Yj ln[l-Pto(0]. (10) 

i Pick / non Pick 

Here the first sum is over all RedMaPPer clusters that have 
been detected by Planck and the second sum is over all RedMaP¬ 
Per clusters that have not been detected by Planck. Figure [T] shows 
the photo-z distribution of the Planck sample (black) compared to 
the subsample (red) defined by the selection algorithm based on de¬ 
tection probability. The data agree in most bins within Icr (of the 
Poissonian errors) and in all bins within 2a. To validate the qual¬ 
ity of the comparison sample we drew 1000 random subsamples of 
250 clusters according to their Pdct and determined the likelihood of 
each subsample. Comparing to the likelihood of the original Planck 
sample, we obtain a p-value of 0.27, so we consider our compari¬ 
son sample as reasonable (i.e., 27% of subsamples have lower like¬ 
lihood than the actual Planck sample). 

We generate a random catalog for the comparison sample by 
using the derived values for the four parameters a^, ft a, a^r and 
ft^ and calculate the detection probability for each entry in the 
RedMaPPer random catalog. 


3.5 Theoretical two-point correlation function 

Our purpose is to compute the cross-correlation between a refer¬ 
ence cluster at given redshift and correlated structures within a cer¬ 
tain redshift range. Note that this differs from computing the usual 
angular correlation function between two samples. In our case, in 
fact, we restrict our calculation of the cluster-group two-point cor¬ 
relation function to redshift bins centered around the reference clus¬ 
ter. The correlated group redshift distribution is thus dependent on 
the reference cluster redshift. We then obtain the total correlation 
function by summing up all the redshift-binned contributions ac¬ 
cording to the cluster redshift distribution. 

The numerical tool we use for calculating the theoretical cor¬ 
relation function is cams souRCESp|Lewis & Challinor|2007|l, which 


* http://camb.info/sources/ 
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computes the angular power spectrum C/S of the matter density per¬ 
turbations, for given input redshift distributions and for different 
cosmological models. We restrict our calculation to standard fiat 
ACDM cosmology {Q.„ = 0.25, h = 0.7) and the linear regime 
only. The relation between the cross-spectra and the projected two- 
point correlation function is given by 



P/(cos 6 ) Cl , 


( 11 ) 


where Pi are the Legendre polynomials of degree /. We use 
a maximum I = 3000 and 8 e [0.01,300] arcmin. We calculate 
the expected two-point correlation (eq. [nj for 20 reference cluster 
redshifts Zci 6 [0.05; 0.5]. For the reference cluster redshift distri¬ 
bution, we assume a Gaussian distribution centered at the cluster 
redshift za, with standard deviation equal to the mean photometric 
redshift error associated to the cluster redshift in the Planck cata¬ 
log, i.e. N(Zci, 0.02). For the correlated groups redshift distribution, 
we use the observed redshift distribution of the RedMaPPer groups 
with richness d > 5, limited to a range of ±0.06, centered around 
Zci. This interval is greater than the bin width in the analysis of 
the observational data of ±0.05 (see section]^, in order to account 
for the errors in photometric redshift (~ 0.02). The observed cor¬ 
relation is the average of the u;, (0) calculated in each redshift bin 
i, weighted by the cluster and group distributions and respective 
average biases, normalized by the total number of objects: 


^(^)lheory “ 


•^20 

2 ji=l 


dA]^ diVf hi fcf u)i{e) 
S^idiV^diVf 


( 12 ) 


Flere &N1 and dAf are the counts per redshift bin of clusters and 
groups respectively. Furthermore, and hf are the average biases 
for clusters and groups within the bin i. We estimate the bias for 
each cluster in the matched Planck catalog and each group in the 
RedMaPPer catalog by using the analytic formula of [Tinker et al.| 
( |2010 ^. This assumes a fixed mass-richness scaling relation, for 
which we employ the result of |Rykofr et al.| ( |2012| l. An analog esti¬ 
mate for the foreground/background structures at |zci - Zgr| > 0.05 
yields a 2pcf consistent with zero within the statistical errors of our 
analysis. 


4 ANALYSIS AND RESULTS 

The relevant quantities we are interested in are the 2pcfs of clusters 
and groups for correlated structure (groups with similar redshifts 
as the cluster), foreground structure (groups with lower redshift 
than the cluster) and background structure (groups with higher red¬ 
shift than the cluster). We compute the angular correlation function, 
where cluster pairs are subjected to one of three constraints: 

(i) the RedMaPPer-Planck cluster pair is separated by less than 
|Az| <0.05. 

(ii) the RedMaPPer-Planck cluster pair is such that Znn < Zpi “ 
0.05. 

(iii) the RedMaPPer-Planck cluster pair is such that Znn > Zpi + 
0.05. 


The first set of pairs allows us to test for the environmental 
impact of physically correlated shuctures, the second for the impact 
of foreground structures and the third for the impact of background 
structures. This selection method is displayed graphically in figure 
with the green volume being the correlated structure, in blue the 
foreground and in red the background. 


We draw 100 sets of 250 Planck random points, assigning 
them the same redshift distribution as the cluster sample. We per¬ 
form the procedure described above on each set, averaging the 
results. We proceed analogously for the comparison sample, by 
drawing 100 samples of unweighted clusters by selecting randomly 
among all clusters of the comparison sample according to their de¬ 
tection probability. The sample size is on average 247, the same 
as the sum over all detection probabilities. The same procedure is 
performed on the random catalog of the comparison sample (see 
section[T4l(. 


4.1 Error estimation 

For estimating errors and covariance matrices, we use three dif¬ 
ferent methods: for the errors of the Planck sample with respect 
to theory, we use a “replace-one” implementation of the Jackknife 
resampling method, for the errors of the comparison sample with 
respect to theory (zero) we use Bootstrap resampling and for the er¬ 
rors of the Planck sample with respect to the comparison sample we 
use a “delete-one” Jackknife resampling by drawing 100 different 
(unweighted) representations out of the complete comparison sam¬ 
ple randomly according to the detection probabilities. The former 
two will be explained in more detail in the following subsection. 


4.1.1 Errors of the Planck sample with respect to theory 

We use a modification of the Jackknife resampling method. In the 
standard ”delete-one“ Jackknife technique, the survey area is sub¬ 
divided into a number of subsamples and the analysis is done a 
number of times equal to the number of subsamples, considering 
each time all samples except one. The Jackknife covariance reads: 

m — 1 ^ 

U ly — ^ (^/,r ~ F/) ~ ’ (13) 

where m is the number of Jackknife samples, is the data 
value in bin i of sample k and x, is the mean value in bin i. Since 
galaxy groups are clustered intrinsically, the errors in neighboring 
bins may be correlated, so we need to take into account the full 
covariance matrix in our analysis. 

We define our Jackknife samples to be equal to the data- 
cylinders in our sub-catalogs. We are using 250 samples containing 
exactly one cluster each and all groups in its vicinity. 

Since our theoretical prediction is made for the exact redshift 
distribution of Planck, we need to find the errors with respect to this 
distribution. A delete-one Jackknife would introduce a systematic 
error here, as the redshift distribution of the sample changes when 
deleting one cluster. To overcome this problem, we use a modified 
Jackknife method: in each Jackknife sample we leave out one sub¬ 
sample (cluster) and assign a weight of two to another cluster. This 
cluster is chosen to be the closest in redshift to the left-out cluster, 
in order to minimize the effect on the redshift distribution of the 
sample. We have to modify equation[T^to account for the changed 
sample size: 

1 

Cij = - V (Xi^k - Xi) (xj,t - Xj). (14) 

The validity of the formula has been verified in a Monte- 
Carlo-simulation. 
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Figure 4. Two-point correlation function for groups in the vicinity of Planck clusters (blue) and groups in the vicinity of clusters in our comparison sample 
(red). In this plot we show the 2pcf for groups with redshift equal to the cluster redshift ±0.05 (correlated structure). The cyan line represents the theoretical 
prediction. Top: complete sample, bottom left: only clusters with S/N >median, bottom right only clusters with S/N <median. In the low S/N case there is a 
slight underdensity in the Planck sample in the region between 10' and 20'. For interpretations see sections [T2]|4.3|and[5] 


4.1.2 Errors of the comparison sample with respect to theory 

In order to perform an error estimation on the comparison sam¬ 
ple with a resampling method, we need a multitude of comparison 
samples. We perform a Bootstrap resampling on the RedMaPPer 
catalog by drawing 1000 random catalogs with the same number 
of clusters as in the original catalog. We then count the number 
of pairs in angular bins around each cluster weighted with the de¬ 
tection probability and compute the covariance in each angular bin 
from these 1000 samples. It turns out that the errors estimated by 
this method tend to be higher than the errors of the Planck sample. 


since we did not account for the modified redshift distribution due 
to the bootstrap here. To overcome this problem we slightly change 
the procedure by bootstrapping sets of 5 groups instead of single 
groups. The sets are created by dividing the catalog into 5 subsam¬ 
ples split by redshift, and selecting one group from each of these 
subsamples. We sort the subsamples by weight, so we ensure that 
each package contains 5 groups with similar weights and different 
(equally distributed) redshifts. In this way the systematic error due 
to the modified redshift distribution is minimized. 
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Figure 5. Same as figurej^ but for groups with redshift Zgr < Zcl ~ 0.05 (foreground structure). The two data sets agree well in the complete sample and the 
high S/N case, while for low S/N a slight overdensity can be observed in the Planck sample nearly over the complete angular region tested, albeit most data 
points still agree within the error margins. 


4.2 Results 

In this subsection we present the results of the angular two-point 
correlation function of galaxy clusters and groups obtained as de¬ 
scribed earlier in this section. We analyze w(6) in 15 equidistant 
angular bins with a width of 4'. We compare the results obtained 
for the Planck sample (blue points in figures [^[^and|^ with those 
for our comparison sample (red points) and with our predictions 
(cyan line). A likelihood analysis is presented in subsection |4.3| 
We expect that a possible effect is stronger for clusters that are 
just above the detection threshold S/N of 4.5. Due to this reason 
we also split the clusters into a high and low S/N sample. The most 


useful approach here would be to split the sample at S/N 7, which 
is the threshold above which the clusters are included in the Planck 
cosmological sample ( [Planck Collaboration et aL||2013dk Unfor¬ 
tunately, in this case the high S/N sample would contain too few 
clusters causing the error limits to become too large, so we decided 
to split the sample at the median S/N 5.4, generating two equally 
large subsamples with 125 clusters each. 

The top of figure|^shows w(ff) for groups at the same redshift 
as the cluster redshift ±0.05 (correlated structure). In the two in¬ 
nermost angular bins both samples are affected by blending effects 
and halo exclusion. The latter is the effect of two nearby structures 
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Figure 6. Same as figurej^ but for groups with redshift Zgr > Zcl + 0.05 (background structure). We observe a slight underdensity in the Planck sample with 
respect to the comparison sample, which is more severe in the low S/N subsample. 


merging into one halo, which has not been included in the theo¬ 
retical prediction. In most bins up to approximately 40' the Planck 
sample shows a slight underdensity with respect to the comparison 
sample, albeit the individual data points still agree within the error 
margins (likelihood analysis shows the underdensity is not signif¬ 
icant, see table [TJ. The excess in the third bin with respect to the 
predicted curve is potentially due to non-linear structure growth. 
In case of the split sample we see a better agreement between the 
two samples for the high S/N subsample (bottom left plot), while 
the agreement is worse in the low S/N case (bottom right) where 
Planck clusters are found in even more underdense background en¬ 
vironments. 


Figure|^shows the 2pcf for groups with redshift Zgr < Zci“0 05 
(foreground structure). The fact that we also observe blending 
here (in the innermost bins), shows that the detection probabil¬ 
ity of RedMaPPer groups also suffers from blending effects, i.e. 
RedMaPPer is less likely to detect groups in the vicinity of a rich 
foreground or background cluster. Besides this effect, one can see 
a slight overdensity in the Planck sample at angular scales >10', 
but the errorbars suggest that this difference is not significant. The 
effect is again weaker in the high S/N and stronger in the low S/N 
subsample. 

Figure|^shows the 2pcf for groups with redshift Zgr > Zci+0 05 
(background structure). Here the 2pcf suffers from blending on 
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small angular scales, too. A slight underdensity can be seen in the 
Planck sample with respect to the comparison sample in nearly all 
angular bins. The individual data points are, however, in agreement 
within the error margins. Also here the observed underdensity ap¬ 
pears less severe in the high S/N and more severe in the low S/N 
case. 

4.3 Likelihood analysis 

In this subsection we investigate the significance by which the 2pcf 
in the Planck sample differ from the comparison sample and from 
the theoretical prediction. We perform a generalized analysis 
that takes into account the full covariance matrix, since as men¬ 
tioned before we expect the errors in neighboring bins to be corre¬ 
lated positively due to the clustering of groups. The generalized 
reads: 

Tgen = • C-ll ■ <5, (15) 

where is the inverse covariance matrix and 5 is the residual 
vector, containing the difference between measured and expected 
values (where measured values correspond to the Planck 2pcf and 
expected values correspond to either the comparison sample or pre¬ 
dicted values) in angular bins. For the foreground and background 
sample we compare the results with zero, since the theoretical pre¬ 
dictions in these cases are several orders of magnitude lower than 
our measurement uncertainty. 

Table[^gives the p-values for all our three different data sam¬ 
ples for Planck with respect to the comparison sample, Planck com¬ 
pared to theory and the comparison sample with respect to the the¬ 
oretical prediction. The four innermost angular bins have been ig¬ 
nored in the calculation, since the data in these bins apparently 
suffer from halo exclusion and blending effects, which have not 
been considered in our theoretical prediction. Thus, the number of 
degrees of freedom is 11. 

The p-values with respect to the comparison sample are typ¬ 
ically quite high (the lowest one being 0.28 for the foreground 
sample), so the null-hypothesis, which states that the two samples 
are similar, cannot be rejected. The p-values are generally slightly 
lower in the low S/N case which supports our assumption that se¬ 
lection effects are predominately observed in the low S/N regime. 
Nevertheless, we cannot confirm a selection bias based on our data 
sample, since the values of the Planck sample and the comparison 
sample are in agreement everywhere. 

When comparing the Planck data with the theoretical predic¬ 
tion, we find high p-values in the correlated and foreground sam¬ 
ples, while we find very low values in the background sample, 
which suggests a selection effect related to lower background den¬ 
sity. To support this assumption, we look at the p-value of the com¬ 
parison sample vs zero (for the background sample) which sug¬ 
gests much better agreement than the value of the Planck sample. 
When looking at the splitted sample with respect to S/N, the p- 
values for the Planck sample are higher than for the complete sam¬ 
ple in both cases, which comes from the larger uncertainties in the 
splitted sample. The p-values for Planck relative to the comparison 
sample in the high S/N case are in good agreement, yet both only 
marginally agree with zero, which we assume comes from cosmic 
variance. In the low S/N sample however, the agreement of Planck 
with zero is significantly worse than for the comparison sample. 
We conclude that the background underdensity for Planck clusters 
is a function of S/N and the effect becomes stronger for low S/N 
detections. 


We also split the group sample in two subsamples at richness 
12 (high A and low A in table[TJ, but we found no significant differ¬ 
ences in these two subsamples. 

Table gives the best-fitting constant values for the 2pcf and 
the corresponding Icr errors. The first four angular bins which suf¬ 
fer from blending have been ignored in this fit. We see that the 
background correlation is not consistent with zero for Planck with 
more than 4cr, while the comparison sample is consistent with zero 
within 1.25cr, which can still be due to statistical fluctuations. Since 
we detect a background underdensity of -0.049 with a significance 
of ~ 4cr with respect to zero but the comparison sample also differs 
from zero with a value of -0.02 at ~ 1.25cr, we conclude that statis¬ 
tical fluctuations in the particular regions used (cosmic variance), 
likely also contribute to the observed defect of Planck background 
groups, but are no sufficient explanation of the full observed effect. 
On the other hand, one could imagine that RedMaPPer detections 
are biased in the vicinity of massive clusters due to the correlated 
structure that surrounds them out to large radii, which might affect 
the detection of groups due to the blending effect, as discussed in 
section l4^ 

When looking at the foreground sample, the slight overdensity 
one might expect from figurej^is not significant, with a p-value of 
0.72. 

4.4 2pcf for Planck and LRGs 

We want to verify our results by comparing them to an independent 
sample of background sources. We replace the RedMaPPer group 
catalog with the CMASS catalog of luminous red galaxies (LRG) 
with spectroscopic redshifts ( [Eisenstein et al.|201 l[|Dawson et al.| 
|2013[[Anderson et al.|2014] l. As clusters and groups tend to feature 
mostly red galaxies, we expect the LRGs to show a similar clus¬ 
tering behavior. Furthermore, if the origin of the underdensity we 
observed is truly the presence of radio sources, which tend to clus¬ 
ter at high redshifts, we expect to see the same effect in CMASS 
galaxies. 

When looking at figure|7] we cannot confirm the underdensity 
that we found for background RedMaPPer groups. If the physical 
effect has a z-dependence, the fact that the redshift distributions of 
the RedMaPPer groups and the CMASS LRGs differ largely might 
be responsible for the observed effect. As we are using uniformly 
distributed random points for the Planck clusters, we do not ac¬ 
count for a potential position dependence of the selection function. 
In particular, a spatial variation of the redshift dependence of the 
Planck detection probability could possibly mimic such a selection 
effect. We investigated the most likely version of this possibility in 
Appendix [a] although more complex dependencies might exist. 

4.5 Implications for SZ and lensing masses 

The potential underdensity in the background of Planck detected 
clusters is nearly constant in all angular bins except the ones that are 
affected by blending. Hence it is straightforward to model the effect 
as the best fitting constant 2pcf in these bins. We use this value of 
-0.049 for further approximations. We estimate the defect SZ sig¬ 
nal caused by this effect in the average beam size of the Planck 
channels that are involved in the cluster detection. We use the me¬ 
dian redshift of the 250 Planck clusters in our sample (0.23) and 
calculate the mean SZ signal of all RedMaPPer groups with red¬ 
shift higher than that value -r0.05 (as our background selection), 
using the scaling relation from|Planck Collaboration et al.H2013d|l. 


© 0000 RAS, MNRAS 000, 000-000 








Environment-based selection effects of Planck clusters 11 


Table 1. P-values for the diiferent samples for Planck with respect to the compaiison sample and to the theoretical prediction and for the comparison with 
respect to the theoretical prediction. 


Sample 


All 

High S/N 

Low S/N 

High A 

Low A 

Correlated 

Plck-Comp 

0.805 

0.555 

0.433 




Plck-Theo 

0.901 






Plck-Comp 

0.28 

0.98 

0.39 



Foreground 

Pick-Zero 

0.72 

0.64 

0.34 




Comp-Zero 

0.34 

0.18 

1.0 




Plck-Comp 

0.48 

0.70 

0.37 

0.89 

0.13 

Background 

Pick-Zero 

0.0060 

0.051 

0.097 

0.18 

0.010 


Comp-Zero 

0.16 

0.023 

0.64 




Table 2. Best fit values and 1 -<t intervals for the foreground and background samples, for Planck and comparison sample. For the background sample, Planck 
is not consistent with zero within 4cr while the comparison sample is consistent with zero within 1.3(T. 


Sample 


All 

High S/N 

Low S/N 

High A 

Low A 

Foreground 

Pick-Zero 

Comp-Zero 

0.040 ± 0.027 
-0.00039 ±0.021 

-0.014 ±0.041 
0.0052 ± 0.033 

0.071 ±0.036 
0.016 ±0.027 



Background 

Pick-Zero 

Comp-Zero 

-0.049 ± 0.012 
-0.02 ±0.016 

-0.047 ± 0.015 
-0.037 ± 0.023 

-0.046 ± 0.016 
-0.017 ± 0.022 

-0.044 ±0.016 

-0.058 ±0.016 


0.15 


0 . 10 - 


0.05 


0.00 


-0.05 


- 0 . 10 , 



20 30 40 
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Figure 7. Same as figure]^ but for CMASS LRGs with redshift zlrg > Zcl + 0.05 (background structure). We do not observe the underdensity that we 
observed for the RedMaPPer group background in figure]^ The small en'orbars come from the much larger number of objects due to the non-overlapping 
redshift distributions of Planck and CMASS (all LRGs are background to all clusters with respect to redshift alone). 
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Analogously we compute the mean SZ signal of the Planck clusters 
themselves. With the number of background groups with respect to 
the previously mentioned Planck median redshift inside the average 
beam size of the involved channels and the average background un¬ 
derdensity we can calculate the fsz -defect caused by a background 
underdensity of -0.049. We find a number of of 1.7 x 10“** rela¬ 
tive to the mean signal of the cluster. This is due to the self-similar 
slope of the Ysz -MOR. Based on this low number, we conclude 
that there will be no implications on cluster masses derived from 
their SZ-signal. For the same reason we conclude the effect of the 
background underdensity on X-ray measurements to be negligible 
as well. 

We estimate the effect of the potential background underden¬ 
sity on the convergence k, which is the quantity that determines 
the magnification in gravitational lensing. As an example, we are 
using a cluster with a mass of M = 3 x 10 '"*Mq, at the Planck 
median redshift 0.23 and sources at redshift of 1.0, assuming that 
5% of the matter between the cluster redshift and the source red¬ 
shift is missing. In a radius of 15’ around such a cluster, the rel¬ 
ative K deficit would then be ~ 80% of the mean k of the clus¬ 
ter. Due to the mass-sheet degeneracy jSchneider & Seitz| |1995[ 
|Seitz & Schneider! 1995[ |1997^ , this large k defect might have a 
much smaller effect on shear measurements. It will, however, have 
a non-negligible effect on magnification measurements. Magnifica¬ 
tion increases the surface area of observed objects with constant 
surface brightness, leading to higher total brightness (lower mag¬ 
nitude). On the one hand, the consequence is a higher (observed) 
galaxy density, as faint background galaxies (just below the detec¬ 
tion limit) might by detected as their brightness increases. On the 
other hand, the increased surface area also increases the separa¬ 
tion between the magnified objects, leading to a lower (observed) 
galaxy density, counteracting the first effect. For steep luminosity 
functions the first effect is stronger (which is generally the case for 
blue galaxies), while for flat luminosity functions the second effect 
dominates (red galaxies). As a consequence, we expect a negative 
2pcf of red background galaxies around clusters caused by this ef¬ 
fect. A potential background underdensity would counteract this 
effect, causing a slightly less negative 2pcf at small angular scales 
and a slightly positive 2pcf at intermediate angular scales. We esti¬ 
mate the amplitude of this effect for a typical Planck cluster using 
equation 10 from |Umetsu et al!] ( |2011| l. We get a result in the order 
of w{6) < 10“^ for 6 = 10' for the effect caused by the magnifica¬ 
tion of the cluster itself, while the counteracting effect caused by 
the underdensity is in the order of w(9) ~ 10“^, both of which is too 
small to be detected in our measurement. 


5 CONCLUSIONS 

Our main scientific goal was to investigate possible selection ef¬ 
fects on SZ selected clusters based on their group environment and 
estimate implications of such an effect on SZ, X-ray and lensing 
mass estimates. 

We summarize our results as follows: 

(i) We do not find an overall selection effect due to correlated or 
foreground structure. 

(ii) We find a potential underdensity of galaxy groups in the 
background of Planck clusters which manifests in an average 2pcf 
in an angular range <40' of -0.049 with a significance of ~ 4cr. 
However, we cannot confirm this effect when replacing RedMaP- 
Per groups with CMASS LRGs in our analysis. 


(iii) This effect grows stronger for low S/N detections and van¬ 
ishes for high S/N detections. We find no dependence of the effect 
on the richness of the groups. 

(iv) We consider three possible explanations for this effect: 

• An erroneous background estimation in overdense environ¬ 
ments might lead to a lower detection probability of low signal 
clusters in these regions. The details and relative importance of 
these effects is likely dependent on the instrumental and survey 
design and the object detection algorithm. On the other hand, the 
fact that Planck detections combine six bands makes this expla¬ 
nation rather unlikely. 

• RLAGNs, which tend to cluster at high redshifts (|Donoso| 
|et al.|2010| [Yates et al.|1989> , contribute to the radio signal in 
regions where the background density is high and could suppress 
low S/N detections. 

• The Planck selection function is responsible for this effect. 

A spatial variation of the Planck selection function that corre¬ 
lates with the spatial variation of the RedMaPPer selection func¬ 
tion could mimic the observed background group underdensity. 
Due to lack of access to the Planck selection function we we are 
not able to test this at this point. On the other hand, we do get 
the same results if we split our sample by distance to the galactic 
disk, as shown in Appendix [A] 

(v) This potential selection effect has a vanishing impact on SZ 
and X-ray mass estimates. The implications on lensing mass esti¬ 
mates are however much larger with an estimated relative k deficit 
of order unity. 

In the latter context, it is interesting to note that|Gruen et al.| 
( |2014^ found a discrepancy from the self-similar slope (yS = 5/3) 
in the i/sz-mass scaling relation for low S/N Planck clusters, with 
a slope of 0.76 ± 0.20. |Sereno et al.| ( |2014| l found a slope of the 
Tsz-mass scaling relation of 1.22 ± 0.24 using all Planck clus¬ 
ters detected by the MMF3 algorithm, and a slope of 1.40 ±0.31 
when using only the cosmological subsample (S/N >7). They made 
an additional analysis forcing the intrinsic scatter to zero, obtain¬ 
ing even lower results for the slopes, 0.95 ± 0.10 for the full and 
1.09 ± 0.17 for the cosmological samples. The background under¬ 
density in Planck clusters that we find in this work potentially ex¬ 
plains their findings, since that could cause a low-biased lensing 
mass, depending on the S/N ratio of the SZ signal, resulting in a 
shallower slope of the scaling relation. The fact that |Sereno et al.| 
( |2014^ find a slightly steeper slope in the cosmological sample, sup¬ 
ports the assumption of the S/N dependence of this effect. [Von der| 
[Linden et al.| ( |2014) , who compared cluster masses from the Planck 
catalog with weak lensing masses from the Weighing the Giants 
project, found evidence for a mass dependence in the calibration 
ratio between the Planck mass Mpi^nct and the weak lensing mass 

M„i which takes the form Mpi^nck ". A possible expla¬ 

nation for their findings might be low-biased weak-lensing masses 
for low-mass clusters, caused by a background underdensity that 
dominates at low S/N, as we hypothesize it in this work. 
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APPENDIX A: SPLITTED SAMPLE WITH RESPECT TO 
GALACTIC DISTANCE 

As mentioned in section [4!4l we found a discrepancy in our results 
as we observe an underdensity in the background of Planck selec¬ 
tions for RedMaPPer groups but not for CMASS LRGs. In order 
to investigate the cause of this difference we splitted the sample at 
the median absolute galactic latitude to find out whether the Planck 
selection function depends on the distance to the galactic disk, as 
it could be caused for example by galactic foreground emission. 
Our uniformly distributed set of Planck random points would not 
account for such an effect. 

The results of the absolute latitude split are shown in figure 
1^ the according p-values are found in table |Al| and the best fitting 
values and l-cr intervals in table |A2| The split results in a nearly 
unchanged result for 2pcf at angular distances up to ~ 30'. Above 
that value however, the underdensity vanishes in the low latitude 
sample while it persists in the high latitude sample. 


APPENDIX B: SPLITTED SAMPLE WITH RESPECT TO 
GROUP REDSHIFT 

As the redshift distributions of RedMaPPer and CMASS are largely 
different, we investigated the possibility of a redshift dependence 
of this effect by splitting the RedMaPPer group sample at z=0.45. 
This value has been chosen to ensure the sample sizes to be ap¬ 
proximately equal for the high z and low z sample after selecting 
the background groups. 

The results of the redshift split are shown in figure [BT| while 
the p-values are shown in table |A1| and the best fitting values and 
l-cr intervals in table |A2[ The split reveals a slightly more distinct 
underdensity of the Planck sample with respect to zero for low red¬ 
shift background groups, as it might be expected from the null re¬ 
sults with the (high redshift) CMASS sample. On the other hand, a 
similar degree of difference between high z and low z background 
can be observed in the Comparison sample, making the relative dif¬ 
ference non-significant. 
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Table Al. Same as table^ this time for the splitted samples with respect to distance to the galactic disk and redshift. 


Sample 

High abs(lat) 

Low abs(lat) 

High z 

Low z 

Pick-Zero 

0.012 

0.25 

0.52 

0.018 

Comp-Zero 

0.17 

0.32 

0.86 

0.0057 


Table A2. Same as table[^ this time for the splitted samples with respect to distance to the galactic disk and redshift. 


Sample 

High abs(lat) 

Low abs(lat) 

High z 

Low z 

Pick-zero 

-0.029 ± 0.023 

-0.023 ±0.016 

-0.037 ± 0.017 

-0.059 ± 0.017 

Comp-zero 

-0.039 ± 0.017 

-0.0023 ± 0.015 

-0.0079 ± 0.0075 

-0.032 ± 0.0099 




Figure Al. 2pcf for background groups with fai* from the galactic disk (left) and close to the disk (right). The sample has been splitted at the median absolute 
galactic latitude. 




Figure Bl. 2pcf for background groups with redshift >0.45 (left) and with redshift < 0.45 (right). 
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